Ferrohydrodynamics: testing a new magnetization equation 
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A new magnetization equation recently derived from irreversible thermodynamics is employed to 
the calculation of an increase of ferrofluid viscosity in a magnetic field. Results of the calculations 
are compared with those obtained on the basis of two well-known magnetization equations. One 
of the two was obtained phenomenologically, another one was derived microscopically from the 
Fokker-Planck equation. It is shown that the new magnetization equation yields a quite satisfactory 
description of magnetiviscosity in the entire region of magnetic field strength and the flow vorticity. 
This equation turns out to be valid - like the microscopically derived equation but unlike the former 
phenomenological equation - even far from equilibrium, and so it should be recommended for further 
applications. 
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The conventional set of ferrohydrodynamic equations 
[1-3] consists of the equation of ferrofluid motion, the 
Maxwell equations, and the magnetization equation. The 
latter was first derived thirty years ago by the author [1] : 



— = f2xM-i(M-M ) 

dt t 
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(Sh'72). Here M stands for the ferrofluid magnetiza- 
tion under the magnetic field H and the flow vorticity 
2f2 = V x v, 77 is the fluid viscosity, <f> — n V the vol- 
ume fraction of magnetic grains in the liquid, n is their 
number density, V the volume of a single particle, and 
r = SnV/ksT is the Brownian time of rotational particle 
diffusion. When the fluid is at rest in a stationary mag- 
netic field, its equilibrium magnetization Mo is described 
by the Langevin formula 

M = nmL(o|, £ = J^, L(Q = coth£ - (2) 

where m is the magnetic moment of a single particle. The 
phenomenological equation (1) generalizes the Debye re- 
laxation equation in case of spinning magnetic grains. 
The spin originates from viscous and magnetic torques 
acting upon the particles - see below Eq. (6). 

Soon after [1], Martsenyuk, Raikher and Shliomis [4] 
(MRSh) proposed another magnetization equation de- 
rived microscopically from the Fokker-Planck equation. 
They have employed for the purpose an original effective 
Held method (EFM) described in detail in [5]. According 
to the method, an arbitrary nonequilibrium magnetiza- 
tion M is considered at any moment as an equilibrium 
one in a certain - specially prepared - effective field H c , 
that is 



M = nmL(()C/(, C = mH c /k B T. 



(3) 



The magnetization (3) relaxes to its equilibrium value 
(2) as the effective field H c (or £) approaches the true 



field H (or £). This relaxation process is governed by the 
equation [4,6] 
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Equations (3) and (4) determine the dependence 
M(f ; H, fi) in an implicit form, dimensionless effective 
field C being the parameter. 

It is well-established that Eq. (4) describes very fine 
real ferrofluids. Comparison with computer simulation 
of the Brownian motion in orientational space [7,8] in- 
dicated that the EFM yields quite accurate results for 
any values of £ and Qt. The same conclusion has been 
made in [6] under comparing of solutions of Eq. (4) with 
the results of numerical integration of the Fokker-Planck 
equation. At that time, the calculations [6-8] have shown 
that the phenomenological equation (1) is valid for any 
field magnitudes but only sufficiently small vorticities, 
fir < 1. The applicability of (1) in the case Or <C 1 was 
corroborated by numerical [9] and analytical [10] solu- 
tions of the Fokker-Planck equation as well. Therefore, 
under consideration of a weakly nonequilibrium situa- 
tion, one should give preference to Eq. (1) as it is far 
simpler for analysis than (4). The latter, however, guar- 
antees the correct description of magnetization even if 
its deviation from equilibrium value (2) is large, Or 3> 1, 
that is when Eq. (1) leads to erroneous results. 

Let us show that a new phenomenological magnetiza- 
tion equation derived recently [11] from irreversible ther- 
modynamics, 
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(Sh'01), is free from the above mentioned shortcoming. 
Namely, being simple enough, the equation is valid even 
far from equilibrium. As a checking on applicability of 



1 



the old (1) and the new (5) phenomenological equations 
we have chosen their predictions about the rotational or 
spin viscosity r/ r of ferrofiuids. Below we shall compare 
i] r obtained from (1) and (5) with its almost exact value 
resulting from the EFM-equation (4). 

A ferrofiuid flow in a magnetic field is accompanied 
with an intertwinement of hydrodynamic and magnetic 
interactions. Being magnetized, the fluid is subject to 
magnetic force and torque with the volume densities 
(MV)H and M x H, respectively. On the other hand, 
the flow vorticity causes a change of magnetization. As 
seen from (f), (4) and (5), fi impedes alignment of M 
with the direction of the local field H. The appearing 
magnetic torque is equilibrated by the viscous torque, 



6n(f>(uj - ft) = MxH, 



(6) 



where u> is a macroscopic (i.e., averaged over physically 
small ferrofiuid volume) angular spin rate of magnetic 
grains. But any deviation of u> from i~2 leads to an addi- 
tional dissipation which is just manifested in rotational 
viscosity. This dissipation contributes to the stress tensor 
[1,2]: 

<Jik = -pSik + V ( J~ + + ^Vp^ikli^l - fij) 

+±.(H i B k -±H 2 5 ik ); (7) 

here eiki stands for antisymmetric unit tensor and the last 
term represents the Maxwell tensor of magnetic stresses 
in ferrofiuids. Eliminating u> — fi from (7) with the aid 
of (6), we obtain [1,2] 
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On recognizing that B k = H k + AixM k , the stress tensor 
(8) takes an evidently symmetric form: 
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as it must be in our approximation. Actually, the ex- 
act result reads as <r ki — o ik — Ie ik i(dui/dt) where 
I = p s d 2 <f>/ 10 is the volume density of the particle mo- 
ment of inertia (p s the particle density and d their diame- 
ter). However, taking into account an extreme smallness 
oil for the particles with d ~ 10 nm, we have not inserted 
the inertia term I(du>/dt) in Eq. (6). 

The boundary wall streamlined by the fluid is acted 
by the force /, = [c^Jrife on unit area; here [ ] denotes 
difference evaluated across the fluid-solid interface and n 
is the normal to the interface. The friction (tangential) 



force exerted on the wall is f T = [cr Tn ] . By using the elec- 
trodynamic boundary conditions [H T ] = and [B n ] = 0, 
we get from (8) 



fr = v(dv T /dx n ) + (M T H n - M n H T )/2 . 



(9) 



For the Poiseuille or Couette flow, v = (0, v(x), 0), in 
a transversal magnetic field, H = (H, 0, 0), magnetiza- 
tion has two components: M = (M x , M y , 0). Then Eq. 
(9) can be written in the form f T = 2(rj + i] r )fl where 
2ft = dv/dx and rotational viscosity is defined as 



7] r = M y H/AQ. . 



(10) 



Thus, the additional viscosity is expressed through the 
off-axis component of magnetization, M y . For small Q,t 
this component is also small: according to the all three 
cited above magnetization equations, M y oc f2r. So, rj r 
does not depend on the flow vorticity in the limit fir <C 1 . 
However, for finite values of fir the viscosity does depend 
on Q. As the result, the function u Tn (fl) deviates from 
the linear one, i.e. a ferrofiuid acquires non-Newtonian 
properties. 

Proceeding to the calculation of rotational viscosity on 
the basis of Eq. (5) , it is convenient to pass from the fields 
H and H c to their non-dimensional values £ and £ : 



dt t ; 2t( 



Cx(Cx£). (11) 



At the stated above arrangement of the applied mag- 
netic field with respect to the fluid flow, Eq. (11) admits 
a steady solution in which the effective field £ tracks the 
true field £ with lag angle a, i.e., £ = (Coosa, ^sina, 0). 
The dependence of £ and a upon £ and fir is given by 
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Substituting M y — nmL(()sina in (10) and using (12), 
we obtain 
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By the same way we find from (4) 
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that results in 
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(13) 



(14) 



(15) 



The solution of Eq. (1) can be presented in the sim- 
ilar form. Let us introduce a new variable £ instead of 
M by the relation M = M (C/0 where M = nmL{^). 
It is worth to note that ( is no more an effective field 
unlike £ in preceding relationships (12)-(15). By substi- 
tuting the components M x — M (C/£)cosa and M y — 
M (C/0 sina in (1), we get cosa = M/M = C/£ and 
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This and the definition (10) yield 
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(16) 



(17) 



Owing to the smallncss of magnetic grains, the Brow- 
nian relaxation time r does not exceed 10~ 4 s even in 
high-viscous fcrrofluids. Hence the inequality fir <C 1 
is usually satisfied. Then, on neglecting the value f2r in 
(12), (14) and (16), all three these relationships are re- 
duced to £ = £. Eliminating now £ from (13) and (17), we 
see that both the old (1) and the new (5) phenomenolog- 
ical equations predict the same dependence of rotational 
viscosity on the magnetic field strength: 
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The EFM magnetization equation (4) yields a somewhat 
different result. Setting ( = £ in (15) gives 
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The viscosities (18) and (19) are compared in Fig. 1. 
Both of them approach the saturation value rj r (oo) = 
3/70/2 at £ » 1. In the figure we plot the reduced ro- 
tational viscosity 77 r (£)/?7 r (oo) as a function of £. The 
upper curve calculated by the EFM [4] represents a 
very good approximation. Actually, as shown in [5,10], 
it hardly differs from the exact solution of the lin- 
earized Fokkcr-Plank equation. Both the phenomeno- 
logical equations, (1) and (5), result in the lower curve 
in Fig. 1 that is described by the Shliomis' formula (18). 
This function agrees with (19) in the low- and high-field 
limits and deviates from that at most on 15 % in the 
entire range of the argument £. 
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FIG. 1. Dependence of the reduced rotational viscosity on 
the dimensionless field strength given by (19) (MRSh) and by 
(18) (Sh'72 and Sh'01). 



When the ferrofluid is subjected to viscous shear, the 
magnetic grains tend to be rotated out of alignment with 
the magnetic field. Thus the flow with a sufficiently large 
shear rate, fir > 1, induces - along with the Brown- 
ian motion - a quotient demagnetization. Formally, this 
effect originates from decreasing the parameter £ deter- 
mined by Eqs. (12), (14) and (16). According to these 
equations, £ = £ when f2r = but the more is fir at 
constant £, the less is £. The reduction of the magne- 
tization leads in turn to some decrease in the rotational 
viscosity. This decrease, imperceptible in practice up to 
f2r ~ 1, becomes then to be very significant. Figure 2 il- 
lustrates the dependence of the viscosity increase on the 
magnetic field strength for three values of the product 
fir. Interestingly, under the finite shear rate the viscosi- 
ties given by (13) and (17) do not coincide with each 
other any more. As seen from the plot, the higher shear 
the more discrepancy between viscosity values predicted 
by the new and the old phenomenological equations. At 
high shear in a high field, Eq. (1) predicts a hysteresis 
of viscosity, which however is corroborated neither by di- 
rect calculations of [6-8] nor by the solution (14)-(15) of 
the EFM-equation (4). Our new equation (5) does also 
not predict such a hysteresis but it provides us with a 
quite satisfactory viscosity description in a wide region 
of parameters £ and fir. Indeed, in this entire region the 
solutions of Eqs. (4) and (5) agree closely, as shown in 
Fig. 2. Thus, Eq. (5) can be recommended for an em- 
ployment on a same level with Eq. (4). It is worth to 
note that all above calculations, carried out for a shear 
flow, apply equally to a rigid rotation of a ferrofluid with 
an angular velocity O in a constant transversal magnetic 
field, H _L f2, and to a quiescent ferrofluid subjected to 
a uniform rotating field H = (H cos fit, iJsinftt, 0) as 
well. 

The difference between Eqs. (1) and (5) also becomes 
manifest if we consider relaxation from an equilibrium 
magnetization in a quiescent ferrofluid after the field is 
suddenly switched off at the moment t = 0. Then the 
fluid remains at rest, fl = 0, so that M and H in (1) 
and (5) are always parallel to H. Hence the equations 
are reduced to 



dM 
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and 
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H e -H 
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respectively. In Fig. 3 we plot the decay of reduced mag- 
netization M(t)/M according to Eqs. (20) with Mo = 
nmL(£) for some initial field magnitudes £ = mff/fcsT. 
The solution of Eq. (20a) reads 



M{t) /M = exp(-t/V) , (Sh '72) 



(21) 



i.e., it does not depend on £. Equation (20b) has the 
solution H c (t) = H- cxp(— i/r), so that we find 



M(i)/M - L(£e-*/ r )/L(0 . (Sh '01) 



(22) 
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FIG. 2. Dependence of the rotational viscosity on the field 
for some values of the shear rate fir, as calculated from the 
EFM [(14)-(15), MRSh], and by the new [(12)-(13), Sh'Ol] 
and old [(16)-(17), Sh'72] phenomenological approaches. 

Thus, only in the limit ( « 1 both decays, (21) and 
(22), are exponential and coincide with each other. Be- 
cause ferrofluids exhibit nonlinear magnetization, the de- 



cay (22) is noncxponential for a finite £ magnitude. Side 
by side with the dependence of viscosity on magnetic 
field and vorticity, this distinct difference in relaxation 
behavior may be of relevance for the interpretation of 
experiments and testing the magnetization equation (5). 




FIG. 3. Time dependence of the reduced magnetization 
M(t)/M after the magnetic field £ is switched off, as de- 
scribed by (22) . The lowest curve represents also the solution 
(21) for any £ values. 
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